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Abstract 

The field theory of nonrelativistic fermions interacting via contact interactions can be used 
to calculate the properties of few-body systems of cold atoms confined in harmonic traps. The 
state-operator correspondence of Non-Relativistic Conformal Field Theory (NRCFT) shows that 
the energy eigenvalues (in oscillator units) of N harmonically trapped fermions can be calculated 
from the scaling dimensions of A^-fermion operators in the NRCFT. They are also in one-to-one 
correspondence with zero-energy, scale-invariant solutions to the A^-body problem in free space. 
We show that these two mappings of the trapped fermion problem to free space problems are 
related by an automorphism of the SL{2, R) algebra of the conformal symmetry of fermions at the 
unitary limit. This automorphism exchanges the internal Hamiltonian of the gas with the trapping 
potential and hence provides a novel method for deriving virial theorems for trapped Fermi gases at 
the unitary limit. We also show that the state-operator correspondence can be applied directly in 
three spatial dimensions by calculating the scaling dimensions of two- and three-fermion operators 
and finding agreement with known exact results for energy levels of two and three trapped fermions 
at the unitary limit. 
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I. INTRODUCTION 



The problem of few-body atomic interactions in the presence of external confining poten- 
tials is motivated by recent advances in experimental atomic physics as well as theory. Recent 



experiments have realized optical lattices with two confined in a potential well 111, 
Such atomic states have been proposed for implementing quantum logic gates [5|, 



Theoretically, the problem of two atoms interactin 
in Ref. Q, see also Refs.0, [TH, H Q T 
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via short-range forces has been solved 
An experimental confirmation of the 
prediction for the ground state energy of two trapped atoms as a function of scattering 
length was recently performed in Ref. [if. If the scattering length of the atoms is tuned to 
infinity and effective range terms are neglected then the S-wave scattering cross section is 
An/p"^ where p is the relative momentum. This cross section is at the upper limit allowed 
by unitarity. The quantum mechanics problem of three particles at the unitary limit in the 
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presence of a harmonic potential has been solved in Ref. [171], see also Ref. 

An outstanding open problem is the many-body problem of fermions at unitary, which has 
been investigated by numerous authors using a wide variety of methods. Gases of trapped 
fermions whose interactions have been tuned to the unitary limit by means of a Feshbach 
resonance have been realized experimentally [TqI . 20, 21|. For a review of experimental and 



theoretical results, see Refs. [2^, [23 



Both the homogeneous unitary Fermi gas as well as 
the unitary Fermi gas in the presence of harmonic traps are clearly of interest. Though the 
many-body physics problem presents physical challenges not present in the two- and three- 
body problems, the existence of exact solutions for N = 2 and 3, where is the number of 

For example, Ref. 
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fermions, can provide important inputs for the case of arbitrary A^ 
proposes a scale-invariant density functional for the unitary Fermi gas whose parameters 
are fixed by matching the known analytic solutions for two fermions at the unitary limit 
in the harmonic trap. Corrections due to a finite scattering length and effective range 
are included in Ref. [25|. This provides another motivation for studying few-body trapped 
fermion problems. 

An interesting theoretical development is the state-operator correspondence which relates 
the problem of finding the energy eigenvalues of A^ trapped fermions at the unitary limit 
to the problem of finding the scalin g di mensions of primary operators in a Non-Relativistic 
Conformal Field Theory (NRCFT) [26|. Note that the NRCFT is defined in the absence of 
an external potential, so the state-operator correspondence relates a property of the theory 
of A^ fermions in free space to the properties of A^ trapped fermions. Another mapping of 



the trapped A^-fermion problem to the free space A^-fermion problem is derived in Ref. [27 
These authors map the problem of harmonically trapped fermions at the unitary limit to 
the problem of finding zero-energy, scale invariant eigenfunctions of the A^-body problem in 
the absence of any external potential. One goal of this paper is to better understand the 
relationship between these two mappings. 

The other main goal of this paper is to show how the state-operator correspondence can 
be applied directly in three dimensions. For two spatial dimensions {d = 2), the theory 
of fermions at the unitary limit is equivalent to noninteracting fermions while in d = 4 
the theory is equivalent to noninteracting bosons j28| . Therefore, in 2 -|- e dimensions and 
4 — e dimensions a perturbation theory in e can be used to analyze the properties of unitary 



fermions [29|, |30 



In Ref. [26|, the e expansion is combined with the state-operator cor- 
respondence to calculate the energy levels of few-body atomic systems in harmonic traps. 
Operator scaling dimensions are calculated in a perturbative series in e and Fade approx- 
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imants are used to interpolate between d = 2 and d = 4 to obtain results for the most 
physically interesting case of d = 3. In this paper, we will illustrate how the state-operator 
correspondence can be applied directly in d = 3. Using a low energy effective field theory 
for two component fermions interacting via S'-wave contact interactions, we calculate the 
scaling dimensions of S'-wave N = 2 and = 3 fermion operators and find agreement with 
the exact solutions for the energy levels of two trapped fermions, as well as the lowest energy 
state of three trapped fermions in an S'-wave. 

The low energy interactions of few-body systems can be studied using the methods of 
effective field theory. These methods are useful when the typical momentum times the range 
of the interactions is much less than one. This is the case for cold atoms, where a complete 
model for the potential is not required and many quantities can be computed in terms of 
the S'-wave scattering length alone. At these energies the S-wave scattering amplitude for 
two fermions with momentum ±p is 

Air 

■ ip 



j[ = 

M pcot 5{p) — 

An 1 
~ M -l/a + rop'^/2 + ...-ip 



where a is the scattering length and rg is the effective range. In the limit r^p <ti 1, effective 
range corrections can be neglected and the two particle scattering amplitude is exactly 
reproduced by a nonrelativistic field theory with a single S'-wave contact interaction. The 
Lagrangian for this nonrelativistic field theory is 

^ = 0^* + ^) ^ " VV^ , (1) 

where is a two-component field operator that annihilates fermion quanta. Here ipip = 
e°''^ipai'i3, SO scattering occurs in the S'-wave, spin-singlet channel only. The coupling con- 
stant, Co(/i), is given by 

where n is the dimensional regularization (DR) parameter and Power Divergence Subtraction 
(PDS) scheme is used to regulate loop integrals j3l|. In DR, loop integrals which are linearly 
divergent when regulated with a cutoff can become finite because DR discards power law 
divergences. These can be restored within the framework of DR by subtracting poles in 
one lower dimension, then the DR parameter, /i, enters the calculation of loop integrals in 
the same way that a hard cutoff would. In the Minimal Subtraction (MS) scheme, where 
the linear divergences are discarded, Co = Aira/M. Thus, it is clear that to make sense 
of the Lagrangian in the limit a ±oo, one needs to use a hard cutoff, PDS, or some 
other regularization scheme that keeps track of linear divergences. We will see below that 
in order to obtain scaling dimension of operators that are consistent with the state-operator 
correspondence of NRCFT, w'e must also use one of these schemes. 

Since the work of Ref. [34|, it is known that the a ±oo limit of the theory of two- 
component fermions in Eq. ([1]) is conformally invariant. ^ The nonrelativistic scale transfor- 



^ For bosons or fermions with more than two degrees of freedom, an S'-wave three-body contact interaction 
is relevant and violates scale invariance 33 1. 
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mation is 



x' = \x t' = \H , t') = A-^/V(^, t) , (3) 

and the nonrelativistic conformal transformation is 

^'-rh '-rh ^•'(^'.'')^(i + c*Pexp(r^^l^),(f,«) . (4, 

It is straightforward to show that the scale and conformal transformations are symmetries of 
the non-interacting theory. For generic values of a the contact interaction in Eq. ([1]) breaks 
these symmetries. Since the two-particle ^-wave cross section is independent of any scale 
when a — »• ±cxo, it is natural to expect the theory to be scale and conformally invariant 
in this limit. Ref. 3^ showed that the off-shell 2 — > 2 scattering amplitude calculated in 



the theory of Eq. ([T]) is invariant under the Ward identities implied by scale and conformal 



transformations when a — > ±oo. Ref. [26j gives a simple argument for why any particle 
number conserving theory that is scale invariant should also be invariant under conformal 
t r ansf ormat ions . 

When a ±oo, the Hamiltonian, H, the generator of nonrelativistic scale transforma- 
tions, D, and the generator of conformal transformations, C, form an SL{2, R) algebra. We 
will show below that the two mappings of the trapped fermion problem to free space fermion 
problems that were discussed earlier are related by an automorphism of the SL{2, R) group. 
This automorphism exchanges the generators C and H. Since the generator C is just the 
external potential for the trapped fermions, this automorphism interchanges the trapping 
potential and the internal Hamiltonian of the gas. Therefore, the automorphism can be used 
to provide a novel group theoretical derivation of virial theorems for trapped Fermi gases at 
the unitary limit. A virial theorem was first derived using the assumption of universality, 
the local density approximation, and thermodynamic arguments in Ref. [s^]. The virial 
theorem was then rederived and generalized using the wavefunctions of the pseudopotential 
model of the unitary Fermi gas in Ref. [23]. Another derivation of the virial theorem using 
the Hellmann-Feynman theorem appears in Ref. j35|. Our derivation is novel in that the 
approach is group theoretical and relies only on the SL{2,R) algebra. The virial theorems 
hold for A^-body energy eigenstates as well for thermal ensembles, and can be applied to 
spin polarized or unpolarized gases. 

The paper is organized as follows: in the next section, we review basic facts about 



NRCFT's and the state-operator correspondence, as well as the correspondence of Ref. [27 
which relates eigenstates of trapped fermions to zero-energy, scale-invariant eigenfunctions 
in free space. We discuss the automorphism of the SL{2,R) algebra which relates these 
mappings and show how it can be used to derive the virial theorems. In section III, we 
derive the scaling dimension of operators with N = 2 and A^ = 3 and show that these 
agree with analytic results for the energies of trapped fermions. In Section IV, we conclude. 
In the Appendix, we solve the problem of two trapped atoms with arbitrary short-range 
interactions. This was first done for arbitrary scattering length in Ref. using the method 
of pseudopotentials. Here we solve the problem by calculating Green's functions for two 
particles in the trap using the field theory of Eq. ([1]). 
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II. NRCFT, 5L(2, R) AUTOMORPHISMS, AND VIRIAL THEOREMS 



In this section we begin by reviewing NRCFT and the two mappings of the problem of 
trapped fermions at the unitary hmit to free space problems j26|, 2^- We then show that 
these two mappings are related by an automorphism of the SL{2, R) conformal symmetry 
algebra. This is the main result of this section. The automorphism is realized by a unitary 
transformation, generated by the Hamiltonian of the trapped fermions, that interchanges 
the internal Hamiltonian of the Fermi gas with the external trapping potential. This unitary 
transformation can then be used to provide a simple derivation of the virial theorems for 
trapped fermions at the unitary limit. 

The many-body Hamiltonian for harmonically trapped fermions in second quantized form 
is the sum of an internal Hamiltonian, ifint, and an external potential, Kxt, which are given 
by 



After the following rescaling. 



V2\ 



2M 



1 



d^x -M uj^ x^ ip^filj . 



(5) 



X 



X 



UJ /i 



which renders all these quantities dimensionless (we are using % = 1 units), we find 



(6) 



int 



ext 



= UJ J d^x -^x"^ ip^ip 
= uC. 



(7) 



Here we have defined Cq^^i) = MColfi), so that Co{^) is independent of M. This shows 
that we can set M = = 1 and measure all energies in units of the fundamental oscillator 
energy, a;. Lengths are measured in units of Oosc = 1/ V Muj. In this section, we will use these 
units and work with H and C rather than i^int and \4xt- C is the generator of conformal 
transformations [26( |. 

If we modify the definition of the scale transformation to include the appropriate trans- 
formation on /i. 



X 



A X 



t' = xH 



A-V ij'{x',t') = X-^/^ij{x,t) 



we find 



H' = X~^H. 



(8) 



(9) 



when a = ±oo. Though the scale transformation of Eq. ([8]) differs from that of Eq. ([3]) by 
additional rescaling of n, the Ward identities derived in Ref. [32] will still hold for any renor- 
malized Green's function that is /j, independent. Likewise, to see the conformal invariance 
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of the theory defined by Eq. ([T]) exphcitly, one must modify the conformal transformation 
in Eq. (jlj) to include a time-dependent rescahng of /i. If D is the generator of scale trans- 
formations then 

H' = e'^^He-'"'^ = e-^^'H {a = log A) , (10) 

which gives the commutation relation [D, H] = 2iH. 

In a NRCFT, the Hamiltonian, H, dilatation operator, D, and conformal generator, C, 
obey the following commutation relations: 

a) [D,H] = 2iH 
h) [D, C] = -2iC 

c) [H,C] = -tD, (11) 

which are the commutation relations of the group SL{2,R). For the theory of Eq. (fl]), we 
have given H and C above and D is given by^ 



'XX ■ j{x) . (12) 



where j{x) is the particle current density. Eq. (fTTb) follows automatically from the defi- 
nitions of D and C and the equal time commutation relations of ip and ip"^. Eq. (fTTb) is 
actually quite general and will hold for any theory in which particle number is locally con- 
served. Note that C = J d^x ^x"^ n{x), where n(x) is the particle density operator. The 
commutator of the Hamiltonian is proportional to the divergence of the particle current [2^ 



[H, n{x)] = —idtn{x) = iV ■ j{x) , (13) 

due to current conservation. Eq. (fTTb) follows by multiplying Eq. ( |T3l) by x^/2 and in- 
tegrating over all space. So if Eq. (ITTk). which is the requirement of scale invariance, is 
satisfied in a particle number conserving theory then the theory will also be invariant under 
the full SL{2, R) conformal group. To complete the algebra of the Schrodinger group (the 
largest space-time symmetry group of free nonrelativistic quantum mechanics), we also need 
the commutation relations of H, D, and C with other symmetry generators: momentum, 
P, angular momentum, J, Galilean boosts, K, and particle number, A^. The nonvanishing 
commutators involving Ki, Pj, D, and N are: 

[K,,Pj]=tN6ij [D,Pi]=tPi [D,Kj] = -tK,. (14) 

Commutation relations involving J are easily deduced from rotational invariance. 
Primary operators in the NRCFT are defined hj O = 0{x = 0,t = 0) and 

[K,O] = [C,O] = 0. (15) 



^ The explicit expressions for C and D are valid at t = 0. For arbitrary t, C{t) = C(0) + t^H — tD{0) and 
D{t) = D{0) — 2tH. The explicit time dependence can be fixed by requiring A{t) — ~i[A, H] + dA/d = 
0,A = C,D, which is required for conserved charges that generate a symmetry of the Hamiltonian. See 
Ref. 



36[ for a one- dimensional conformally invariant quantum mechanical system. 
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Note that primary operators are defined to be located at the origin of space and time. The 
particle number (No) and scaling dimension (Ao) of the primary operator are defined by 

[D,0] = lAoO, 

[N,0] = NoO. (16) 
If we translate the primary operator O to another point in space-time 

then it is straightforward to show using the commutation relations listed above that 

[K,0] = {-itd, + NoXi)0 

= -i{fdt^tx-d^tAo)0^%^NoO. (18) 

The field has = — 1 and = where d is the dimensionality of space. The density 
operator, ip'^ip, has N^]^ = and A^t,/, = d. For a finite conformal transformation we have 

0'{x,t) = e-'^^0{x,t)e'^^ 

^^PhT7T^^K^^^>^^ • (19) 



{i + xt)^o ^V2(i + At)y \i + xt'i + xt 

which agrees with Eq. (j4]) for the case O = ip. 

Next we discuss consequences following from the algebra in Eq. ffTTj) . Let us define 



-f^osc = H + C 

L± = H-C±iD. (20) 

The L± are ladder operators that raise and lower energy eigenvalues of i^osc by two oscillator 
units, as can be seen from the commutation relations 

[L±,Hosc\ = =F2L± 
[L_,L+] = 4i7o3c, (21) 

which are easily derived from Eq. (fTTl) . Eigenstates of ifosc come in infinite towers of equally 
spaced states. The ground state of one of these towers is denoted by IV'o) which satisfies 
L-\^o) = 0. 

The problem of finding energy eigenstates for the trapped particles can be mapped to 
the free space theory in one of two ways. The method of Ref. |26|] begins by noting that 

e^L_e-^ = -C, (22) 

and furthermore 

CCt|o) = [C,O^]\0) 

= 0, (23) 

where is a primary operator and |0) is the vacuum. Then 

L.e-^O^O) = -e-^CO^O) = (24) 
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so the ground state of the tower is lipo) = e ^O^O). We find that from any primary operator 
we can construct a tower of eigenstates of i^osc- It is straightforward to show that [26| 

HoscC-^O^O) = e-^{C - iD)O^\0) = Aoe-^O^|0) . (25) 

Thus the scahng dimension of the operator in the NRCFT (in the absence of an external 
potential) gives the ground state energy (in oscillator units) of the corresponding state \ipo) 
in the problem with an external harmonic potential. 

Ref. [23| pointed out another correspondence between eigenstates of the trapped fermions, 
i^osc, and zero-energy, scale invariant eigenstates of H. The result of Ref. [13] can be obtained 
starting with a relation analogous to Eq. 



e^L_e-^=H. (26) 

From this relation it is clear that \iPq) = e~'"\ipj), where is a zero-energy eigenstate 
of the Hamiltonian, H\ipi,) = 0. For this state to be an eigenstate of Hose, it must be an 
eigenstate of iD as well, 

tDm = m- (27) 



The energy eigenvalue of \ipo) is 



= e-"" {H + iD)\^lj,) 
3 



z/+-iV)|^o). (28) 



To understand the significance of z/, note that the A^-body wavefunction associated with the 
state \ipi,) is 



TV 



i'M) = (oin^(^"'*)IV'.)- (29) 



i=l 



Using 



it is easily seen that 



e-*"^V^(xl)e'"^ = ei"7/;(e"xi) (30) 



i'.yjj =A-'MSi), (31) 

so the A^-body wavefunction for the state is a homogeneous function of the A^-body 
coordinates. Note that the iV-body wavefunction for the trapped problem is given by 

{xi\ipo) = {xi\e~^\ip^) 

= e-S«^"/2^,(f,). (32) 
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To understand the relationship between the two mappings of the trapped problem to free 
space problems, we observe that Eqs. (122|) and (126!) are related by the following automor- 
phism of the SL{2, R) algebra 



H ^C. 



D 



-D. 



(33) 



This is an automorphism of the SL{2, R) algebra which is implemented by a similarity 
transformation using the elements 

Jn{n+l/2)Hosc 



9n = e 

whose action on the generators of SL{2, R) is ^ 









f c 


9n 1 


C 




H 











(34) 



(36) 



These identities immediately lead to the virial theorems for trapped fermions at the unitary 
limit derived in Refs. 27|, |3J]. The thermal expectation value of an arbitrary operator, O, 
is given by 

(O) = Tr[e-^(^°--'^+^+-^-^-)0] , (37) 

where we have included separate chemical potentials, ^+ and yU_, for spin up and spin down 
fermions, respectively, so our results can be applied to spin polarized as well as unpolarized 
gases. For the expectation value in Eq. ([ST]), or the expectation O in an eigenstate of i/osc 
we have (gnOdn^) = {O), because [Hosc,N±] = 0. Therefore, 

(if") = (C"), (D2"+1) = 0. (38) 

For n = 1, this implies (-f^osc) = Eq = {H + C) = 2(C) which is the virial theorem first 
derived in Ref. [s^]- The generalization to arbitrary moments of C in the ground state is 
straightforward: 



71—1 



L:+L. 



I^o) 



2 
Eo 
2 

Eq 



(V-olC^-^l^o) + ^(^o|2(n - 1) C^'^ + [H, C 
(V-olC^-Vo) + 7(^o|2(n - 1) C"-i + C 



Eo + in- 1) 



4 



(39) 



Eq. ([36]) is a special case of 









/ 


cos^e* 






{'] 








cos^ 9 








v 


sin 29 


- sin 29 




(35) 



The automorphism in Eq. 



is obtained for sinl 



±1. 
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We have used L-\'iI)q) = = {iI)q\Lj^. This simple recursion relation immediately gives all 
higher moments of the trapping potential which can be written in closed form as [27| 

= . (40) 

This concludes our general discussion of NRCFT. The main result of this section is the 
automorphism of SL{2, R) which relates the two known mappings of the trapped A^-fermion 
problem to problems involving the N fermions in free space. This automorphism provides 
a simple, group theoretical method for deriving virial theorems for both eigenstates and for 
thermal expectation values with arbitrary chemical potential for the two spin components. 
In the next section of the paper, we will show that the state-operator correspondence can 
be used directly in d = 3 using the effective field theory of Eq. ([1]). 



III. STATE-OPERATOR CORRESPONDENCE IN d = 3 
A. Two fermions 

The problem of two fermions interacting via short-range interactions in the presence of an 
external harmonic potential is exactly solvable ||^ . This solution is reviewed in the Appendix. 
The ground state energy of two fermions at the unitary limit in a harmonic trap is 2, in 
oscillator units. In this section, we verify the state-operator correspondence by evaluating 
the anomalous dimension of the composite operator ipip using the NRCFT of Ref. [32|. We 
compute the matrix element 

{Q\Z^4fi)^Pm-p), (41) 

which is given by the Feynman diagrams in Fig. [H The factor Z^^{^) is required for 
composite operator renormalization. We work in the center of mass frame, where E is the 
total kinetic energy and the momentum of each particle is p = \p\ = V ME. The diagrams 
form a geometric series 

{o\zUfi)^^P\p,-p) = .42) 



where 

G^(0,0) 



{2tiYe-P/M 



M 



-^i/"- V-p' + ^e) . (43) 



The first line of Eq. (143!) is obtained after evaluating by contour integration the energy 
integral in the one-loop bubble graph that is pictured in Fig. [H Note that the one-loop 
graph in the NRCFT is related to Green's function for the free two-body Hamiltonian: 



E — Hq 



2) J {2tiYe-P/M 



(44) 
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FIG. 1: Feynman diagrams contributing to the matrix element in Eq. (|41|) . The black blob is the 
operator ipip. 

where d is the number of spatial dimensions and the factor (/x/2)^~'^ is inserted to give the 



correct dimensions. We use DR and the PDS scheme SlJ to evaluate the integral. The 
integral is linearly dependent on n in the PDS scheme, reflecting the linear divergence, but 
IJ, independent in the MS scheme. Keeping the linear /j, dependence is critical for finding the 
correct anomalous dimension for the composite operator ipip. 
The result for the matrix element is then 



(/i) 



l + fCo(/i))(/i + *p) 



4:71 f Z^^{fx)\ . ^^^^ 



M \ Coifi) J 1/a + zp' 

where we have used Eq. (j2]). The matrix element is /i independent if Z^^{fi) oc Co(yu) and we 
can fix the constant of proportionality by demanding {0\Z^^{fi)iljtp\'p, —p) = 1 for = 0."^ 
Then 

ZM^^) - 

= (46) 
1 — /i a 

The anomalous dimension of the the operator ipip is then given by 

= . (47) 

1 — /xa 

The effective field theory is a NRCFT when we take the limit a —>■ ±oo, and then 7^^ = — 1. 
The scaling dimension of ipip is the naive dimension, 2A^ = 3, plus the anomalous dimension, 
7i/)i/) = —1 so A^,^, = 2A^ + 7^^ = 2, in agreement with the state-operator correspondence. 
Note that one must take the logarithmic derivative with respect to fi at finite a, then take 
the limit a —>■ ±00. If the limit is taken prior to computing the derivative, then Z^^{fi) = 0. 
However, this is an artifact of the boundary condition that {0\Z^^{fi)iljilj\p, —p) = 1, which 
is no longer possible when a = ±00. If we start with Eq. (H5|) . take the limit a — ±cx) 
we should demand that residue of the 1/p pole be a /i-independent constant and we again 
obtain 7^^ = — 1. 



^ The same normaUzation condition is obtained if one requires that the sum of ah Feynman diagrams yields 
the same resuh when evaluated in the minimal subtraction (MS) scheme, in which case all loop graphs 
are finite and Co = Ana/M. 
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Another way of obtaining the scahng dimension of a primary operator, O, in NRCFT is 
to consider the two-point function: 

Go{x,t) = (O|0(f,t) 0^(0,0)10) 

oc e{t)t-'^° exp (^-iNoY^ , (48) 



where the second hne of Eq. (1481) follows from scale and Galilean invariance 2^ . Note we 
assume Nq > 0. Fourier transforming this Green's function yields 



Go{p,E) = J d'^xdte-'^'^'P-'^Goix,t) 



oc ^ (49) 

This formulae can be used once additive renormalizations are carried out. For example, in 
the noninteracting theory, for d = 3, 



G^4p, E) = -i^ - ^J-ME + ^ + ^e] (50) 

After removing the yU dependence using an additive renormalization, we can compare with 
Eq. (H9|) to obtain A^^ = 3, which is the correct answer for a free theory. For the interacting 
theory, 

Goip,E) = I ^-— + ^ = 1 . (51) 

In deriving this result we have included the factor Z^^{fi) computed earlier. The first term 
can be removed by additive renormalization or else we can remove the cutoff dependence by 
taking fi ^ oo. Then, the second term then yields A^^ = 2 when comparing with Eq. fH^ 
in the limit a — > ±cxd. 

An alternative formulation of the NRCFT employs a composite field, which in the context 
of nuclear physics is called the dibaryon formalism [s^]. This formalism is used most often 
in three-body calculations. In the dibaryon formalism one introduces a composite field, (j) 
that has the same quantum numbers as iIjiIj and removes the four-fermion interaction using 
a Hubbard- Stratonovich transformation |37|. In this formalisrn, Eg. (H9l) can be directly 
compared to the dibaryon propagator (see, e.g., Eq. (3) of Ref. |37|), in the hmit a — >• ±cx), 
ro — and again one finds A^ = A^^ = 2 

It is also interesting to see how the scaling behavior of the two-body wavefunction dictates 
the anomalous dimension of the corresponding two-body operator in the field theory. This 
sheds further light on the relationship between the results of Ref. [23] and Ref. [2y]. The 
unrenormalized sum of all graphs in Fig. [1] can be expressed quantum mechanically as 

(oi^^Ip, -p) = J-0yi {q\i + ^^J\P) ' (52) 

where T is the transition operator that is a solution to the Lippmann-Schwinger equation, 

T = V + V^^T, (53) 
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and matrix elements of V are {q\V\p) = Cq. Here, we have let the number of spatial 
dimensions, d, be arbitrary. It well known from nonrelativistic quantum mechanics that the 
exact solution to the scattering wave equation with incoming particles with momentum p is 

X,ix) = {x\l + ^^J\p) (54) 

so we can write 

{0\i/;m-p) = XpiO), (55) 

so the matrix element can be interpreted as two-body wavefunction at the origin. However, 
this is divergent for the interacting theory. For two fermions at the unitary limit the two- 
body wavefunction, Xp{^)y is proportional to r^^*^ for small r. Inserting a complete set of 
state into Eq. (!52|) and regulating the expression with a hard cutoff in momentum space, we 
obtain 

{0\^l^m -p) = xp<0) = J''^ I d'^xe-^'^-'xpix) . (56) 

This integral is of course divergent. The degree of divergence is determined by the x — 
behavior of Xp{^) which is independent of p. It is easily seen that the integral diverges like 

If renormalize the matrix element in Eq. (!52l) with a multiplicative factor of Z^^(A), we must 
have Z^^(A) oc h?^'^ to get a finite answer for the matrix element. This leads to 7^^ = 2 — d 
which gives for the scaling dimension for A^^ = d + 7^^ = 2, which is the correct answer 
for arbitrary d. 



B. Three particles 

The three-body problem in the presence of an external harmonic potential with infinite 



two-body scattering length was solved in Ref. |17| . For any interaction to take place two of 
the three particles must be in an S'-wave. Ref. 17|] solved the three-body problem for arbi- 
trary /, where / is the total angular momentum of the three-body system, using the method 
of pseudopotentials. Since the interaction is modeled as zero-range, the three particles are 
free except when the coordinates of two of the particles coincide. The wavefunction is then 
a solution to the free Schrodinger equation subject to the boundary condition (for arbitrary 
a) 

lim V'(fi, f2, f'i) (X — -- + 0{rij) (sj = -Sj) (58) 

where Vij = \fi — rj\, Si and Sj are the spin quantum numbers of particle i and j, respectively, 
and the limit Vij ^ is taken holding the coordinate of the third particle fixed. (For Si = sj 
the wavefunction must vanish as Vij — »■ 0.) In this paper, we will only consider the case of 
1 = 0. It would be interesting to extend the analysis to arbitrary I but that is beyond the 
scope of this paper. 
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First we briefly review the solution obtained in Ref. 17|. Suppose we choose the spin 
states so that si = S3 = —S2. The three-body wavefunction is parameterized as 

^(fi, f2, fs) = (1 - P13) V^cm(^cm) ^ p) , (59) 

where Rem is the center-of-mass coordinate, r = |ri — r2|, p = |2r3 — ri — f2\/V^, and the 
operator P13 interchanges ri and f^. The wavefunction for the center-of-mass coordinate, 
"00111 (-Rem); is a solution of the simple harmonic oscillator Hamiltonian, so i?cm = uji'in + Z + 
3/2). The function xi^^ p) obeys the following differential equation 

+ ^ - ^(-^ + P^) + MiE - Bern)) Xir, p) = 0. (60) 

Imposing the boundary conditions in Eq. fl58l) . and demanding the wavefunction be finite as 
p 0, one finds that 

|:x(o, p) + Ixio, p) - [^P. j = , X(r, 0) = . (61) 

For a = ±00, it is possible to solve the boundary condition with a factorized solution, 
x(r, p) = F„(_R)0„(a), where 2 R'^ = + p^, and a = arctan(r/p). The function 0„(a) is 
determined by 



i2 



-0„(a) = So„0n(a) 



0„(7r/2) = 0, 

0^(0) = -^M^m ■ (62) 



while F{R) obeys the differential equation 

^2 T a r,2 



In addition we must have 0„(O) 7^ 0, so that the residue of the l/r^ pole in Eq. ( l58l) is not 
equal to zero. The first two lines of Eq. ( l62l) are solved by 



a oc sm 



TT 



a - - ) So,n 



(64) 

while the third line of Eq. fl62p leads to the transcendental equation for so,n 



so,n COS I — — ) + ^ sm I — — 1=0. (65) 



V 6 

Note the solutions come in pairs, So,n = ±|so,n|- So,n = ±2 is a solution to Eq. (!65!l . however, 
inspection of Eq. (!64l) shows that for sq,™ = ±2 , 0„(O) = 0, which will not satisfy the 
boundary condition in Eq. (l58l) . There are no other integer solutions to Eq. (165|) . and all 



14 



remaining solutions to Eq. (165|) give nontrivial solutions to the three body-problem. Numer- 
ical values of the five smallest values of \so,n\ are 2.16622, 5.12735, 7.11448, 8.83225, 11.06273. 
The numbers so,n determine the energy eigenvalues via Eq. (!63l) . The solutions of Eq. (163|) 



are 
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F„{R) oc R'"'" e"^'*^"/' L['^'-\R^Mu) , (66) 

where Lg^°'"'* is a generalized Laguerre polynomial, and the energy eigenvalue is E = + 
'^(so,n + 1 + 2g). Note that for the wavefunction to be square integrable, we must have So,n 
positive in Eq. fIBB]) . The dependence on the quantum number q shows that for each so,n 
there is an infinite tower of evenly spaced states whose energies are separated by 2uj, as 
expected from the SL{2, R) algebra. 

At this point we would like to demonstrate the correspondence between trapped eigen- 
states and zero-energy, scale-invariant eigenfunctions of the free Hamiltonian. To find these 
states, we can choose the same variables, -Rem, R, and a, which were used to solve the 
trapped three-body problem. The function ^cm(-Rcm) is now a solution to the free particle 
Schrodinger equation, ipcm{Rcm) oc e"*^'^™'^™, and we should take Pcm = to obtain a zero- 
energy state. The eigenvalue equations for are still Eq. fl6Tl) . and Fn{R) obeys Eq. (|63|) 
with uj = E = Ecm = 0. Thus, the solution for the zero-energy scale-invariant wavefunction 
has Fn{R) oc R^°'" and the zero-energy, scale- invariant solution to the three-body equation 
is 

Mn, rs, rs) oc (1 - P^) — 

r p 

= (1-Pi3)i?^°--^^^. (67) 
sm(2a) 

Clearly, the scaling exponent for this state isu = so,n — 2, so Eq. fl28l) tells us that the energy 
of the ground state of the infinite tower of states is E = uj{5/2 + s^n), in agreement with 
the result obtained by direct solution of the three-body equations. 

Now we would like to see how the effective field theory reproduces these results. We will 
show that the effective field theory allows one to derive a bound state equation which exhibits 
scaling solutions whose scaling exponents yield energy eigenvalues via the correspondence 



of Ref. [271]. Then we study how the state-operator correspondence can be used directly in 
three dimensions by calculating the anomalous dimension of an operator in the NRCFT and 
verifying that it reproduces the known result for the lowest energy state of three trapped 
particles in an S'-wave. 

In applications of effective field theory to three-body problems it has been found useful to 



employ the dibaryon formalism discussed earlier [33] . In the present context the composite 
field should be called a difermion, which we will denote 0, which has the same quantum num- 
bers as i/jip. A Hubbard-Stratonovich transformation is used to trade the contact interaction 
in Eq. ([T]) for (p^ipip and ilj^ip^(f) couplings. Loops of ipip contributing to the d) self-energies 



are summed to all orders to obtain the propagator. We refer readers to Ref. [33|] for details 
on this procedure. 

The scattering of (p and proceeds via an infinite number of ladder-like diagrams. These 
can be resummed using a one-dimensional integral equation, which is pictured in Fig. [2l 
Double lines are (p propagators and single lines are tp propagators. Evaluating the diagrams 
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+ 




FIG. 2: Integral equation for (pip scattering. 



in Fig. [2] and projecting onto the S-wave yields the half-off-shell integral equation 



S{p,k) 



M 
~2pk 
1 

71 



log 



p'^ + pk + k'^ - ME 
p"^ -pk + k'^ - ME 



(68) 



7 1 , ( p^ +pq + 
dq — log 



pq 



ME\ 



pq + q^-MEJ -1/ a + ^?,q^ / A - ME 



Siq,k). 



Here the momentum k is on-shell, 3A;^/4 = ME + and the momenta p and q are off- 
shell. Note that S{p, k) corresponds to the sum of ladder diagrams and does not include the 
LSZ factors required to obtain the on-shell amplitude when p = k, nor is it normalized to 
give the three-body scattering length for p = k = 0. Apart from these factors, the equation 
obtained here for (pip scattering is identical to that obtained in Ref. |33i] for three nucleons 
in the J = 3/2 state of nucleons. Ref. 33|] defines a half-off-shell amplitude, a{p), which 
is normalized to the three-body scattering length, 03 = a{p = k). The function S{p,k) is 
related to the function a{p) of Ref. [33| by 



S{p, k) 



3M 



a{p) 



8 l/a+ ^^p^/A- ME 



(69) 



It is straightforward to reproduce the results of Ref. [17l| using the NRCFT and the 



mapping of Ref. [27|]. To find a zero-energy, scale-invariant eigenstate of the free space 
problem we can consider the equation for three-body bound states in the limit a —>■ ±00 and 
E = 0. The bound state equation is obtained from Eq. (!68l) by dropping the inhomogeneous 
term in the integral equation. In the limit E and a — ±00, the bound state equation 
becomes 



S{p, 0) 



TT 



dq 
p 



log 



p"^ + pq + q^ 
p2 _ pq _^ q2 



S{q,0) 



(70) 



Then one looks for solutions of the form S{p, 0) oc p~^o,n-i^ which is possible if So,n satisfies 
Eq. (165!) (33|. The integral equation for S{p, k) is finite and does not require renormalization. 
In order for the second diagram on the right hand side of the integral equation in Fig. [2] to 
converge for large g, S{q, 0) must vanish as g — 00, which then leads to S{p, 0) oc 
To see how the p ^ 00 behavior of S{p, 0) is related to the scaling behavior of the many- 
body wavefunction, we recall the three-bo dy p osition space wavefunctions can be obtained 
from S'(p, k) using the following transform [38| 



x{f, p) 



p sm \\l\pp 

dpS{p, k)- 
-l/a+Jlp^-ME 



-r-\J\p^-ME 



(71) 
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FIG. 3: Diagrams contributing to the renormalization of the operators Oi. 



It is straightforward to show that the function x{i"^ p) obeys Eqs. ( |60f6Tll . Eq. ( |60l) and the 
second boundary condition in Eq. fl6T]) follow directly from the definition in Eq. flTTl) . while 
the first boundary condition in Eq. fl6Tl) can be obtained using the integral equation for 
S{p,k). Taking the limit E = 0, a = ±00, and inserting the the asymptotic solution for 
S{p, 0), we find 



x{r, p) oc I dpp ^ sin ^ ^ ^p R cos a | e~ 



3/2 pR sin a 



OC i?"*" " sin 



TC 

S0,n I " - 2 



(72) 



which is the correct form of the zero-energy, scale-invariant solution. 

Finally we wish to understand the state-operator correspondence for the case of the three 
trapped fermions in an S-wave. As an example, we show how the state-operator corre- 
spondence can be used to calculate the lowest energy state of three harmonically trapped 
fermions in an S'-wave. We compute the scaling dimension of the operator 



(73) 



Operators that are a total time or space derivative are not primary, so Oi is the unique 
primary operator with one time derivative. S'-wave operators with two space-derivatives can 
be put in the form (pV'^ip after integration by parts, and are therefore equivalent by the 
equations of motion for ip. Therefore, in the noninter acting theory, where has dimension 
3, Oi is the unique operator with naive dimension 13/2. Note that the lowest energy state of 
three noninteracting fermions in an S'-wave has energy 13/2 u, so the naive scaling dimension 
is consistent with the state-operator correspondence for the free theory. We compute the 
matrix element {0\Zi{A) Oi\p, —p) . The diagrams that contribute to this matrix element 
are pictured in Fig. [3], which shows a tree-level graph and another graph which includes 
the half-off-shell amplitude S{p,k). The off-shell legs of S{p,k) are contracted with the 
operator Oi to form a loop. This graph sums all loop corrections to the matrix element 
(0|Zi(A) Oi\p, —p). The sum of all diagrams contributing to the renormalization of Oi in 
the limit i?, 1/a = is given by 



Zi{A) 



—2 



2po^S{p, 0) 



dpp^ S{p, 0) . 



(74) 



We have included the factor Zi{A) for composite operator renormalization. A cutoff on the 
virtual loop momentum is used to regulate the loop integral. The asymptotic form of S{p, 0) 
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is determined by Eq. (ITOll . so in general we have 



^(p,0) = 5^c„p-l^°-l-^ (75) 

m 

where the coefficients in the expansion, c^, must be determined numerically from the solution 
to the full integral equation. Obviously the large p behavior is dominated by the smallest 
values of m in Eq. (!75l) . For the operator Oi, the only divergent contribution comes from 
the term m = 1, |so,i| = 2.16622. All other terms in Eq. (!75l) give a UV finite contribution 
to Eq. (17^ . The integral on the right hand side of Eq. (174|) diverges as A'^"'**"-^', so 

7i = aAzi(A) = |so,i|-3. (76) 

The scaling dimensions of and ip are = 2 and = 3/2, as discussed in the previous 
section, and the time derivatives add 2 to the naive dimension of Oi. Adding the anomalous 
dimension, 7^^, we find the scaling dimension, Ai = 5/2 + |so,i|, which, via the state- 
operator correspondence, is also in agreement with the result for the lowest energy state of 



three harmonically trapped fermions in an S-wave [17| . 

We should not consider the operator (pip. The analog of this operator in the formulation 
of the theory without a difermion field would be {ipip)ip = e'^^ipaippip which is not allowed 
because of Fermi statistics. A local operator which creates (or annihilates) three fermions 
at a point must have derivatives acting on at least one of the fermion fields. Therefore 
the operator (pip which seems allowed if (p is treated as a boson, must be excluded from 
consideration when classifying local operators in the NRCFT. This can also be seen from the 
state-operator correspondence for the noninteracting theory. The naive scaling dimension of 
(pip is 9/2 in this case, but there is no state of three trapped fermions with energy 9/2 uj. So 
clearly one obtains a contradiction with the state-operator correspondence if (pip is allowed. 

Note that the true ground state of three fermions at the unitary limit has / = 1. It 
would be interesting to derive the transcendental equations analogous to Eq. (!65l) for / 7^ 
from the integral equations for scattering in higher partial waves derived in Ref . [3^] . These 
should give the numbers si^n that determine the energy eigenvalues of three fermions in 



higher partial waves [17|. Another problem is to determine operators that correspond to 
states with energy eigenvalues E = 5/2 + |so,n|, n > 2. These come from operators with 
two or more time derivatives or four or more spatial derivatives, or mixed time and space 
derivatives. In the equations analogous to Eq. (1711) . these operators will lead to more factors 
of po or which will make the integral more divergent. This leads to more terms in the sum 
in Eq. fl75l) contributing to the anomalous dimension. It should be possible to find a basis 
of operators in which the anomalous dimensions are given by the so,n that are solutions to 
Eq. ([65]). 



IV. CONCLUSIONS 

We have studied the problem of fermions with infinite two-body scattering length con- 
fined in harmonic traps. The A^-body problem can be mapped to problems involving A^ 
fermions in the absence of an external potential. One approach is to map solutions of the 
trapp ed p roblem to zero-energy, scale-invariant solutions to the Schrodinger equation in free 
space [l7| . Another approach is to relate the energy levels of A^- fermion states to the scaling 
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dimensions of primary operators in an NRCFT [26*]. In this paper, we have shown that 
these two mappings are related by an automorphism of the SL{2, R) conformal algebra of 
the NRCFT. This automorphism interchanges the internal Hamiltonian of the NRCFT with 
the harmonic trapping potential. This provides a simple, group theoretical way of deriving 
virial theorems for trapped Fermi gases at the unitary limit. The virial theorems apply 
for energy eigenstates as well as thermal ensembles and hold for both spin polarized and 
unpolarized gases. 

One goal of this paper was to apply the state-operator correspondence [26] directly in 



the three spatial dimensions {d = 3), which is clearly the most important case. In Ref. [26 
the state-operator correspondence was combined with e expansions about d = 2 and d = 4 
to do perturbative calculations of the energy levels. We sought to apply the state-operator 
correspondence directly in three dimensions using the NRCFT of Eq. ([1]). This is clearly 
more difficult because analytic results are only available for two fermions. For two fermions 
we showed how to use the state-operator correspondence to calculate the energy levels of 
two trapped fermions at the unitary limit. For three fermions, the NRCFT gives an integral 
equation for (pip scattering which can used to find the zero-energy, scale-invariant eigenfunc- 
tions which can be used to find the eigenf unctions of the three trapped fermions via the 
correspondence of Ref. [13] • We showed how to use the state-operator correspondence to 
calculate the energy of the lowest energy S-wave three-fermion trapped state. It would be 
interesting to extend application of the state-operator correspondence to all eigenstates of 
the trapped three fermion problem. 

Since the problems of two and three trapped fermions in the unitarity limit can be solved 
using quantum mechanics and the pseudopotential boundary conditions of Eq. (155]) . an 
important question is whether the mappings of the trapped fermion problems to free space 
problems will be useful for obtaining new results. The virial theorems jl^] are an example of 
results that the conformal symmetry of the NRCFT can provide in the absence of an exact 
solution of the quantum mechanics problem. It would be interesting if the integral equations 
of the effective field theory could be used to calculate corrections to energy levels from a 
finite scattering length, or obtain new results for problems with four or more fermions at 
the unitary limit. It would also be interesting if SL{2, R) invariance can be used to obtain 
information about correlation functions of two-point functions of primary operators in the 
eigenstates of harmonically trapped fermions. For example, if SL{2, R) invariance provides 
interesting constraints on correlation functions like Eqs. fl48ll49p with the vacuum replaced 
by the ground state of N trapped fermions, one could perhaps learn something about the 
low lying excitations of the ground state of a trapped gas of cold atoms at the unitary limit. 
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APPENDIX A: TWO FERMION IS A HARMONIC TRAP IN EFT 



In this appendix we solve the problem of two-particles interacting via short range forces in 
the presence of a harmonic potential. This problem was first solved in Ref. m and is typi cally 
analyzed using the method of pseudopotentials, see e.g. Refs. 0, Q Ef H H 13, 0- 
Here we solve it by evaluating the two-particle Green's function. 

Consider the Green's function 

'^i?tot(^3,^4;^l,^2) = (X3,f4| ^^ ^ff(0) l^l'^2), (Al) 

where H{H^) corresponds to the interacting (noninteracting) Hamiltonian. 
GEtoti^Sj ^2) obeys the integral equation 

^Stot(^35 3^4; a^i, 2:2) = G%^^^{x3,X4;xi,X2) 

+ Coif^) j d^yG%^^{x3,X4;y,y)GE,,,{y,y;^i,^2) ■ (A2) 

This equation can be derived in quantum mechanics using a delta-function potential with 
coefficient Gq{^) or from the Feynman diagrams of the field theory in Eq. ([T]) in position 
space. It is helpful to go to center of mass coordinates 

Xi,2 = /?±V f3,4 = i?'±^r', (A3) 

because the Hamiltonian factorizes in these coordinates. The noninteracting Green's func- 
tion is given by 

(^3, 0:4, XI, X2) - E,ot - E{n) - E{m) ' ^^^^ 

Here n = {n^, ny,n.^) and rh = {rrix, my, rriz)- The iprn{R) are eigenfunctions of the Hqm with 
energy E[rh), and (f)^{R) are eigenfunctions of the noninteracting ifrei with energy E{m). 
The interacting Green's function G^^^^ has the same form as G^tot with (l)%{R) replaced by 
eigenfunctions of ifrei, (pn{R)- Since Hem is the same in either case, the iprniR) are common 
to G^^^^^{xs, X4; Xi, X2) and G£;^^^(x3, X4; Xi, X2), we can project onto an energy eigenstate of 
Hqm- Defining 

Gg)(r',f) = J rf'^^t/'^^>^(^)V'^(^')G'gi(^3,X4;fi,f2), (A5) 

where E = i?tot ~ E{m) and E{m) is an eigenvalue of i/cmj we find that Ge{t' ., r) obeys the 
equation 

G^(f', f ) = Gl{r\ r) + CM G%{r', 0) Ge{0, r ) . (A6) 

The notation is similar to that used in Eq. fj44l) . however here Ho is the simple harmonic 
oscillator Hamiltonian. Setting r ' = we obtain 

G^(0, f ) = ^"^^"'"^ ^ ^ . (A7) 

l-Co(/i)GUO,0) 
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The poles of this expression are solutions to 



— -GUm-O, (A8) 
where (7^(0,0) is the Green's function for the simple harmonic oscillator, and is given by 

POO 

G0,(0,0) = - / rft(0|e(^-^o)*|0) 
Jo 

\ Etf Mu , 



A transcendental equation similar to Eq. (lASP but with a different regulator was obtained 
in Ref. 40]. This integral is evaluated for negative E — Eq, where Eq is the ground state of 
the oscillator, using dimensional regularization: 



Jo {smhty/^ Jq ^ 



,_,r[i-f]r[f + f] 



= 2^^/2-1^ 2il2 4J _ (AlO) 

r[i-f + f] 

Analytically continuing the integral from negative to positive E we find, for arbitrary d, 

(All) 

Just like the G^(0, ) in the absence of the oscillator potential, this integral is linear diver- 
gent, but finite if evaluated using dimensional regularization. The integral is defined exactly 
as in the free space theory, multiplying the integral by {^/2Y~'^ and subtracting the pole at 
c/ = 2. We find 

G^(0,0) = -(-/i+v^^if-^j . (A12) 
Therefore we find the poles of the Green's function are located at. 



3 _E ] 

2uj\ 



4vr I a V\\- £] 



(A13) 



which is the transcendental equation first derived in Ref. . This result is easily generalized 
to include effective range corrections. Effective range corrections and higher order terms in 
the effective range expansion can be incorporated using higher dimension operators with 
derivatives. We can choose a basis where each higher dimension operator contributes a 
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factor of C2n{ME)"' to the tree level scattering amplitude, see Ref. [4l| for more details. 
The scattering amplitude is in the absence of an external potential is 



A 



(En^2„(^)(ME)")-^ + f(^ + Zp) 

Att 1 



Mp cot 6{E) - ip 



(A14) 



Including the higher derivative operators in the Eq. flA6l) for the Green's function, one finds 
the the formula in Eq. flA13l) becomes 

= -(-pcot^(E)-v^ ^|t_Sj I ' (A15) 

In Ref. [l3 it was pointed out that Eq. (IA13I) receives significant corrections when aylMuj = 
o/oosc > 1- Later, it was shown [ill. Il2i. ilS] showed that reliable results could be obtained 
by making the substitution 

^^-p cot 5(E). (A16) 

This substitution was called the "effective-scattering length model" , which we see here can 
be derived in a straightforward way using effective field theory. 
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